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Abstract 

The possible existence of Higgs-Higgs bound states in the Higgs sector of the 
Standard Model is explored using the \hh) + \hhh) variational ansatz of Di 
Leo and Darewych. The resulting integral equations can be decoupled exactly, 
yielding a one-dimensional integral equation, solved numerically. We thereby 
avoid the extra approximations employed by Di Leo and Darewych, and we 
find a qualitatively different mass renormalization. Within the conventional 
scenario, where a not-too-large cutoff is invoked to avoid "triviality," we find, 
as usual, an upperbound on the Higgs mass. Bound-state solutions are only 
found in the very strong coupling regime, but at the same time a relatively 
small physical mass is required as a consequence of renormalization. 

PACS numbers: 11. 10. St, 14.80.Bn, 11.80.Fv, ll.10.Gh 
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I. INTRODUCTION 



The possible existence of bound states for the Higgs boson has been studied by several 
authors |l]-|J with both perturbative and non-perturbative calculations. At present, there 
is little agreement between the quantitative predictions of such calculations. 

A variational method, within the Hamiltonian formalism ||, has been used by Di Leo 
and Darewych (DLD) Q. However, because of the apparent complexity of the resulting 
integral equations, they resorted to additional approximations that are unsatisfactory when 
the coupling is strong. In this paper we show that an exact decoupling of DLD's integral 
equations is possible by use of some symmetry properties. This allows a considerable sim- 
plification of the problem. For an s-wave solution we show in detail that the method gives 
rise to a one-dimensional integral equation that can be tackled numerically. 

Mass renormalization (beyond normal ordering) plays a crucial role, although it is fi- 
nite. (The infinite mass renormalization in DLD turns out to be an artifact of their other 
approximations). The physical mass is significantly reduced from its classical value in the 
strong-coupling regime. Because of this mass-reduction effect we find that the occurrence of 
bound states (i.e., solutions in the 2-particle sector with E < 2m) is shifted towards the very 
strong coupling regime, well beyond the reach of any perturbative approximation. However 
a relatively small physical mass (m m 0.1 — 0.5 TeV) is required as a consequence of the 
same renormalization effect. 

Our paper is organized as follows: the variational method is described in section II. A 
general prescription for the mass renormalization is then provided and discussed. Section III 
deals with the delicate aspect of mass renormalization through a variational one-particle trial 
state calculation analogous to the two-particle trial state used in section II. The existence 
of Higgs-Higgs bound states is discussed in section IV and the results are compared with 
those of other authors. 

Our discussion here uses the conventional framework that the Higgs theory is an effective 
theory with a large, but finite cutoff A J7|. 



II. VARIATIONAL METHOD 

Following DLD, our starting point is the Lagrangian 

£ = -id^h - \m\h 2 - ^Xvh 3 - ^\h\ (1) 

for a neutral, scalar Higgs field h. For the present we regard mo, v, and A as three indepen- 
dent bare parameters; only in section IV shall we need to impose the constraint ra\ = |At> 2 
that arises when this model is obtained from a spontaneously broken A$ 4 theory. 

In the Schrodinger representation, at t — 0, the field is quantized in terms of creation 
and annihilation operators 

/i(x) = / = akexp(zk • x) + aj c exp(— zk • x) (2) 



/ (27r) 3 2w k 

satisfying the usual commutation relations 
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= 6 3 (k-p). (3) 

The energy of the single particle states is 

u(k) =u k = Vk 2 + m 2 (4) 

where m is the physical mass which may differ from the classical mass mo. The Hamiltonian 
H is obtained from Eq. (|l|) as 

d 3 k \ uj k - ^— J a k a k - J d 3 k ^— (a k al k + a k a_ k ) + 

^ jd z xh 3 + jjd z xh^ 



+ : 



(5) 



where 5m 2 = m 2 — m\. The Hamiltonian has been normal ordered with respect to the 
physical mass m. 

The preceding discussion follows the conventions of Ref. ||], except that our A is a factor 
of 6 larger than theirs, since we use A/4! rather than A/4. 
The trial two-boson bound state considered by DLD M is 

|tf 2 ) = Jd 3 pB(p)al > al p \0) + Jd 3 pd 3 qd 3 kG(p,q,k)alal l ai\0)6 3 (p + q + k), (6) 

with |0) the vacuum annihilated by a. We observe that the function B(p) may be taken to 
be symmetric without any loss of generality; a general B(p) could always be decomposed 
into an even and an odd part and the odd part would give no contribution to This 
reflects the fact that a bound state of two identical bosons must be even under spatial 
inversion. Similarly, we must have G(p, q, k) = G{— p, — q, — k). Furthermore, there is no 
loss of generality in assuming that G(p, k, q) is invariant under any permutation of the three 
momenta p, k and q. (Of course, because of the momentum-conserving delta function, the 
function G really involves only two independent momentum arguments.) 
The B and G functions are determined from the variational principle 

5(ty 2 \H-E\ty 2 } = (7) 

which provides two coupled eigenvalue equations: 

5(^ 2 \H - E\^ 2 ) 



6B*(k) 

5(^ 2 \H - E\^ 2 ) 
5G*(p,q,k) 



(8) 



(9) 



Explicitly, these equations are: 

^ m2 \ A r 3 B(p 



' om z \ . . A f B(p) 

2u k E Bk + — — d 3 p^^ + 

V u k J QAix^uJk J u p 

3A?; r r 1 

+ tw / d3 p / ^ ^ 3 ( k + p + q)) , g ( k > p. q) = ( 10 ) 
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3 / Sm 2 \ ^ 
A 3 



G(k 1 ,k 2 ,k 3 ) + 



Xv 



(4!)7r 3 / 2 v /^iu;2W3 ~i 

1/2 



647T J ^^2^3 ~i 



£ MJ d 3 q 5 3 (k t + p + q) f -^-) G(k u p, q) = 



(11) 



where uJi = u k .. These equations are equivalent to Eqs. (8), (9) of Ref. |||, but their structure 
appears considerably simpler because we have taken advantage of the symmetry properties 
mentioned above. We can achieve an exact decoupling of these equations by the following 
algebraic manipulations. First, we introduce the auxiliary function 



A{k) =u k J d 3 p J d 3 q 5 3 (k + p + q 



G(k,p,q 



(12) 



so that the first equation becomes 
3\v 



-A(k) = [E - 2u k 



5m 2 



B , ] s x f d3 P 
\ ) 647r 3 u; fe - 



/^B(P) 

J UJrt 



(13) 



Similarly, re- writing its last term in terms of A(kj), Eq. ([ll]) becomes 
G(k u k 2 M) 



-A 1 Ef =1 


'A{ki) + (87r 3 / 2 /3)^(k 4 )" 


t>47T 3 y/UJx UJ2OJ3 




— E 





(14) 



Note that the resulting form of G manifestly respects the symmetry properties invoked 
earlier. Inserting this equation back into Eq. (|i~2"|), the second eigenvalue equation Eq. ( |i~T|) 
is equivalent to 

A(k) + (8tt 3 / 2 /3) i J ^}^ vB(k) + 



+ 



2lui 



1 + J(k) 

d 3 p [ d 3 qd 3 (k + p + q)K(k, p, q) \A(p) + (8tt 3/2 /3)vB(p 



1 + J(k) 

where the kernel K(k, p, q) is wholly symmetric: 

A 



(15) 



^(ka,k 2) k 3 



64vr 3 



UJ1UJ2U3 



(16) 



and 



J(k)=cu fc Jd 3 pK(k,p,-k-p) 



(17) 



is a logarithmically divergent integral. We shall regularize it with an energy cutoff, 
Vp 2 + rn' 2 < A. 



The eigenvalue equations (fTq) , ([151) may now be easily decoupled by replacing A from 
Eq. flT3p into Eq. (13). In this way, we obtain the following integral equation for i? 
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E-2u 




J(k) 



ri/i \ A 1 + 3 J(k) /-d 3 p 



1 + J(k)) 



l + J(k), 

J d 3 p K(k,p, -k - p)B(p) (5m 2 + \v 2 + uj p (E - 2uj p )) . (18) 



For an s-wave 5 function, the angular integration can be performed analytically, yielding 
a one-dimensional integral equation for B which can be solved by numerical methods. There 
are two main conceptual problems that must first be dealt with (i) the mass renormalization 
parameter 5m 2 needs to be determined, and (ii) the integral J is logarithmically divergent 
and requires regularization, say with an energy cut-off A. This last point has to do with 
the physical interpretation of the whole theory. The current orthodox viewpoint is that the 
original A$ 4 theory is only an effective theory, valid up to some finite energy scale A that 
acts as a cutoff. A is then another parameter of the theory, in addition to mo and A. We 
shall adopt this viewpoint in this paper. (For a heterodox viewpoint, see Ref. 0.) 

Mass renormalization is crucial since the existence of bound states hinges on the com- 
parison between the energy E and the energy of two free bosons at rest, 2m. Any attractive 
self-interaction that tends to bind two particles will also give rise to a reduction of the phys- 
ical free-particle mass compared to the classical mass m . Thus it would not be legitimate 
to ignore mass renormalization and just impose m = m . 

The form of the left-hand side of equation (|18|) suggests @] that the desirable mass 
renormalization is such that the combination 



should vanish. Since m 2 should not be k dependent, we define 



5m 2 = -\v 2 J{ °] (20) 
1 + J(0) V ; 

For an infinite cutoff, J — > oo and we would get 

m 2 = ml - Xv 2 (21) 

which is a finite mass renormalization. [In DLD, due to their other approximations, the 1 + 
J(0) denominator is absent in Eq. ( pOf) , so that they found an infinite mass renormalization.] 
In the next section, we will show that the above mass renormalization is justified by 
considering a variational calculation of a one-particle state. 



III. SINGLE-PARTICLE MASS RENORMALIZATION 



The mass renormalization prescription (|20| ) may be recovered in an analogous self- 
consistent variational procedure. In this case, the trial state j^k) for a single boson with 
mass m and momentum k is taken to be 

|* k > = C(k)4|0) + |rf 3 pD(k,p)ai +k a t _ p |0) (22) 
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where D(k, p) = -D(k, — p — k). The variational principle now requires that 

5(y k \H - E(k)\V k ) = (23) 
which leads to the coupled eigenvalue equations 



°w(*-£-*«)+i5./ 



Xv f d 3 pD(k, p) 







(24) 



£>(k,q) 



5m 2 ( 1 



E(k) 



+ 



A y d 3 pD(k,p) 

647T 3 7 ^^q^k+q^p^p+k 

At; C(k) 

167T 3 / 2 ^k^q^k+q 



+ 



0. (25) 



Eliminating the integral we obtain 
D(k,q) 1 



^k 



C(k) 



i )7r 3/2 v 



^q^k+q 



^q + ^k+q - ^ 



1 



£(k) 



Substituting back in Eq. 



then for we find 



where 



Jo(k) 



A 



64tt 3 



£(k) = ^ k 



1 



2^k 



8m + A?/ 



■/o(k) 
1 + J (k) 



^p^k+p 



^p + ^k+p — ^T" 



f ^ + - J -) - £(k) 



(26) 



(27) 



(28) 



This has a similar structure and the same ultraviolet behaviour as the integral J(k) of Eq. 
(pi). In principle, the presence of the single particle energy i£(k) in the denominator requires 
us to solve Eq. (|27j) and Eq. ( p8|) self-consistently. However, when A ^> m (as it should be) 



J (k) « J(k) 



A 



A 



ln(A/m) + finite terms 



128vr 3 7 |k| 3 32tt 2 
In such a limit, we self-consistently obtain E(k) = provided we take 

Jo 



(29) 



5m = —\v* 



l + Jo 



(30) 



to be compared to Eq. (p0|) . Once more for A>mwe recover the mass renormalization 
prescription (pT|). 

Full consistency would also require that when E(k) = uj^ the function D should vanish. 
In fact, since is the energy of a single particle state |k) = a k |0) with mass m, the state 
l^k) can have a mass m and a Lorentz-covariant dispersion relation only when |\I>k) = |k). 
This requirement is not entirely trivial: the mass renormalization prescription (|30"D is just 
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what we need in order to guarantee that the single particle state |k) effectively is the lower 
energy one-particle state for the full interacting Hamiltonian (g). Inserting Eq. (|27| ) and 
(p0|) into Eq. (p6|), we find that the ratio D/C — > for A — > oo as we expected. In other 
words that means we find a quantum-field renormalization constant Z = 1. 

It is instructive examining the same result from the point of view of standard perturbation 
theory. If there were some "rule" forbidding the existence of states with more than two 
particles in the Fock space, then the variational trial state ( ]2"2]) would lead to an exact 
result. The same result should be achievable by perturbation theory provided that we sum 
all Feynman diagrams whose intermediate states do not contain more than two particles. 
However, in the self-consistent variational procedure the mass of the single particle state |k) 
is supposed to be the true physical mass m ^ m . In the language of perturbation theory 
this is equivalent to associate a renormalized propagator with each internal line of Feynman 
diagrams. The easiest way to do that is by re-writing the Lagrangian (U) as 

C = C + C 1 + C 2 (31) 

where Co is the zeroth-order non-interacting part 

Cv = -\d^h R d»h R -\m 2 h 2 R (32) 

C\ is an interaction part contributing at tree level 

d = -\{Z - l^hnd^hn - \{Z - l)m 2 h 2 R + \Z5m 2 h 2 R (33) 

and £2 is the interaction 

C 2 = ~XvZ^h\ - \z 2 \h R (34) 



Here h R is a renormalized field h R = h/yZ. Imposing that the renormalized propagator 
has a pole at p 2 = — m 2 with unit residue gives two conditions || : 

Z5m 2 = -HL p (-m 2 ) (35) 

and 

z = i + ^nL> 2 )U 2 (36) 

where p 2 = p 2 — u 2 , and i(27i) 4 U* oop (p 2 ) is the sum of all one-particle-irreducible diagrams 
containing loops. Such diagrams can only arise from the lagrangian part C 2 - In our reduced 
Fock space the only diagrams contributing are the bubble-chain diagrams reported in Fig.l. 
These are all naively divergent but, if a regularization prescription allows their resummation, 
their infinite sum yields a finite contribution 



where the one-loop term is 



11 r 1 n2 



-- 2 z 2 xi{p 2 ) 



-^Z 2 XI(p 2 ) 



+ ... (37) 
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iJU„> 2 ) = ~(^W)/(p 2 ^ 



and I(p 2 ) is the divergent integral 



Hp 2 



d 4 q 



1 



i(27r) 4 (q 2 + m 2 — ie) ((p + g) 2 + m 2 — ie) 
The infinite series (^) can be exactly summed up 



(38) 



(39) 



nr oo > 2 ) = (z\v 2 ) 



^Z 2 I(p 2 ) 
1 + \\Z 2 I{p 2 ) 



(40) 



In order to mantain Lorentz-covariance, the integral ( p9|) can be evaluated by dimensional 
regularization. By use of the Feynman formula and Wick rotation, in a d dimensional 
Euclidean space the integral reads 



Hp 2 



71 



d/2 



(2vr) 4 

For d = 4 + e and e — > we obtain 
1 1 



r(2-d/2) / m 2 +p 2 x(l -x 



8?r 2 e 16vr 2 
Insertion in Eq. (HUf) gives 
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+ hi7r + y In (m 2 + p 2 x(l — 



x 



+ 0(e) 



^n* oop (p 2 )|_ m2 ~e 2 

and from Eq. ([36|) then 

Z = 1 + C(e 2 ) 

while from Eq. (|35|) the mass renormalization parameter reads 

A/(167r 2 e) 



5m = —\v' 



A/(167r 2 e) - 1 



(41) 



(42) 



(43) 



(44) 



(45) 



Thus, in the physical d = 4 space, we recover for e — > the mass renormalization prescription 
of Eq. (H) and Z = 1. 



IV. SEARCH FOR HIGGS-HIGGS BOUND STATES 

In the previous sections we described the variational method and discussed its internal 
consistency from a general point of view. So far the three parameters ra 2 ,, v and A have been 
viewed as independent and we have not specialized to any particular physical problem. We 
now wish to use the method described in section II to search for Higgs-Higgs bound states 
in the Higgs sector of the electroweak theory. 
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The scalar sector of the electroweak theory has the form: 

£ = - A($2 _ rff + const . ( 46 ) 

Defining the 'Higgs' field h by $ = v + h, we obtain the original Lagrangian ([!]) with the 
parameters related by 

ml = hv 2 . (47) 

It can be shown that the same relation holds for v being the minimum of the Gaussian 
effective potential 0. In the standard interpretation, we also have a large but finite cutoff 
A. The theory is approximately Lorentz invariant for energies small compared to a finite 
energy cut-off A. The vacuum value v is fixed empirically in terms of the Fermi constant 
G F 

2v 2 = (48) 

Thus the bare mass m Q is proportional to the square root of the coupling A. The physical 
model is entirely described by two independent energy scales: the energy cut-off A and the 
bare mass m which also fixes the coupling strength through Eq. (fl7|). 

For a large cut-off A, according to Eq. ([?!]), we expect a mass correction 5m 2 « —Xv 2 



which overcomes the tree-level mass (f|7]). Since m 2 is positive definite we expect that 
something should prevent it from becoming negative. In fact the general discussion of section 
II must be modified when the physical mass m approaches zero. At the point m — > 0, A — > oo 
the integral J(0) is not analytical, and some extra care is required in handling the two limits. 
As previously discussed J(0) diverges logarithmically according to Eq. fl29"|) for any finite m, 
while it vanishes linearly for m — > at any fixed cut-off A. Thus for any large but finite A the 



coupled equations ( |T7|) and (pOj) must be solved together yielding a real cut-off dependent 
mass m(A). Of course at any finite cut-off such equations mantain a fc-dependence since the 
theory is not Lorentz invariant. We define m as the k = value corresponding to the energy 
required to create a boson at rest. For k = a generic scattering solution of the integral 
equation ( jTB| ) has E = 2m where m is determined from Eq. (P0| ) by insertion of Eq. fl4"7|) 



2 2 

m = m r 



2.7(0) 



1 + J(0) 

while J(0) is the integral (|i"7D of the kernel (|l^) evaluated at ki = and at E — 2m: 



(49) 



J (0) = ^-J 9 J x (50) 

32tt z Ji or — ax — p 

where g = A/m, a = (3 — m\j 'm 2 ) / '4 and (3 = (1 — m^/m 2 ) /2. 

The two coupled equations ( f49D and (|50|) give the physical mass m. The integral in Eq. 
([501) can be evaluated analytically: 

3tt? 2 
62n z v z 
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where 



f(g)=hi(g+y/g*-l) + -'£XJ^=-]n(l-g 7± + Jg*-ly/>£-l 



- \ E ln (i - 9i± - , (52) 

and 7± = (a ± v 7 ^) /2 and A = a 2 + 4/3. 

A numerical solution for m as a function of the coupling parameter m is reported in 
Fig. 2 for several cut-off values. In the weak-coupling limit, for m smaller than 0.3 TeV, 
we recover the perturbative solution m ~ m which holds up to a quite huge A. That 
is equivalent to neglecting J(0) altogether in Eq. ([49]). For larger couplings the solution 
deviates from the perturbative regime and the physical mass m is heavily reduced and 
strongly dependent on the cut-off choice. Notably in the range 1 TeV < mo < 2 TeV, 
where several authors find Higgs-Higgs bound states, the physical mass is spread over a 
large energy range going from m = mo for A = mo to m rs mo/100 for A = 35 TeV. However 
the upper bound m = mo is only reached for an unphysical cut-off equal to the mass, which 
makes the integral J vanish. A cross-over is observed at a critical A c = 3.05 TeV from 
a monotonic increasing behaviour of m versus m , to a non-monotonic beahviour with m 
rising to a maximum and then decreasing for larger m . The maximum value of m never 
exceeds the critical value m c = 1.1 TeV (which is reached at a coupling m = 4 TeV) for 
any choice of coupling and cut-off. Thus, for any A > 3.05 TeV and any m , we always find 
a physical mass m < 1.1 TeV. 

Moreover, in the strong coupling regime the physical mass becomes very small: in such 
strong-coupling limit a simple analytical solution may be found by requiring that m <C m 
and expanding both equations (|49|) and (|50|) in powers of m 2 /m§. By use of the analytical 
expression (p2"|) , the integral (0) may be written as 

J(0) = + °^ 2 ^ ( 53 ) 



which to first order in m does not depend on the coupling mo- Eq. ( [49]) may be inverted 
and to the same order in m 2 reads 

J(0) = \ + 0(m 2 /m 2 ) (54) 

Eliminating J(0) yields 

4 t> 2 

which is consistent with our assumption that m mo provided that mo ^> 47r 2 t> 2 /(3A) or 
explicitly m • A 3> 0.8 TeV 2 . When such condition is satisfied the physical mass does not 
depend on the coupling and tends to a finite limit proportional to A -1 . This behaviour is 
evident in Fig. 2. 

For completeness we should mention that whenever mo > A a second larger solution for 
the physical mass m comes out from the coupled equations (^9|), (fjCf), but such solutions 
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probably have no physical meaning. We remark that the range mo > A could be of physical 
interest since the cut-off A should be compared to the physical mass m which is generally 
significantly smaller than the coupling parameter m . 

Returning to the bound-state problem, let us insert the numerical solution of the coupled 
equations (pE9|), fl5"T|) into the integral equation (|T8|). Neglecting the slight k-dependence of 
J(k) m J(0) (which is the only approximation we are making apart from the choice of the 
trial state) we obtain the following integral equation 



(2u k - E)B(k) = - J d 3 p}C(k, p, -k - p)B(p), 



(56) 



where the kernel /C is defined as 

1 / 2ml + m2 * 



/C(k,p,q) 



'5ml — 2m 2 ' 
2m 2 , + m 2 



1 



m 2 + 2ml + u p (E — 2uj p ) 



UJ k + Up + UJ q + 



J_ 



E 



(57) 



We notice that the second term inside the brackets contains an attractive part plus an energy 
dependent part proportional to (E — 2uj p ) which is always repulsive if E < 2m, and thus 
weakens the bonding of any bound state. It is instructive to see how the approximate integral 
equation of DLD @ can be recovered from our almost exact treatment of the variational 
method. That can be done by taking the perturbative limit (J(0) = 0, m = mo) and by 
neglecting the energy dependent repulsive part proportional to (E — 2u p ) in the second term 
of the kernel according to their approximation E w 2uo p . Moreover, since they also take 
E Up + uj k oj ( 



+ Uk, the denominator of the second term is approximated as 



1 



0J k + U p + 0Jq 



E 



1 



yielding the approximation 



1 



UJ„ 




(58) 



/C(k,p,q) 



A 3m 2 

64-7T 3 UJ k 0Jp 



3m 2 



U>* UJqUJp 



(59) 



which is Eq. (16) of Ref. 0. This approximation requires m ps mo, which is valid only in 
the region mo < 0.3 TeV, according to our numerical results in Fig. 2. DLD found bound- 
state solutions with the kernel (|59|), but only for larger couplings mo > 0.9 TeV; this is 
well beyond the perturbative regime HJ9|JT0[| , and well beyond the region of validity of their 
approximations. 

To address the problem of the existence of bound-states beyond the perturbative regime, 
we must use the full integral equation (56|). The most interesting case is an s-wave solution, 
and in that case the integration over angles may be carried out exactly yielding a one- 
dimensional integral equation. We change the integration variables according to 



d 3 p... 



2tt 



Ut 



ujpdujp 



UJqduJg . . . , 



(60) 
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where 



+ uo 2 — m 2 ± 2 



\/{ul-m 2 )(uj$-m 2 ). 



Then integrating over uj q gives 



,A 

(2u k -E)B(k) = J r (u k ,u p )B(p)du p , 

J m 



(61) 



(62) 



with 



2m\ + 777,2 ( / S 777 ^ — 2m 2 



lQir 2 v 2 uJk 1 \ 2m 2 , + m 2 



/Wp — m 2 + 



2ml + m 2 + u p (E — 2u p ) 



777/ 



In 



LU i — fL 



UJ- 



in 



(63) 



where f2± are the poles of the kernel K in the variable uo q 



E — u>k 



m 



+ 



1 

± ~\ 
2\ 



E -uj k 



UJr, 




2m\ + 2m 2 . 



(64) 



The one-dimensional integral equation ( |B"i2"D may be solved numerically by standard matrix 
techniques. (One must first choose the parameters mo and A, and find the corresponding 
physical mass m from the coupled equations (|49|) , (|5T|) .) We can describe three different 
scenarios: 

i) A < A c = 3.05 TeV (small cut-off). Since m is a monotonic increasing function of the 
coupling strength mo, both m and mo must be small compared to A. In this weak-coupling 
limit there are no bound-state solution and the lower eigenvalue of Eq . (|62D is the free particle 
energy E = 2m. 

ii) A > A c , m < 2 TeV (large cut-off and moderately strong coupling). Beyond the 
critical value A c = 3.05 TeV the Higgs mass m is not a monotonic increasing function of the 
coupling strength mo (see Fig. 2). This strength can be large since m is bounded and never 
comparable to A. Beyond the weak-coupling limit, where there are no bound-state solutions, 
an intermediate range can be described for m w 1 — 2 TeV. In this range non-perturbative 
effects are evident since m is a decreasing function of the bare mass mo (as shown in Fig.2). 
In this range bound state solutions have been found by several authors ]I]-f|]. However the 
strong reduction of the physical mass m, in comparison with the bare mass mo, rules out 
the existence of bound states with E < 2m in this regime. 

iii) A > A c , mo > 2 TeV (very strong coupling). For very large mo the renormalized Higgs 
mass m saturates at the finite value given by Eq.([55|) (see also Fig.2). A further increase 
of the coupling strength allows the occurrence of bound-state solutions whose precise onset 
depends on the chosen energy cut-off A. In Fig. 3 we show the binding energy (E — 2m), 
in units of 2m, for a typical cut-off A = 4 TeV, just above A c . An eigenvalue E smaller 
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than 2m appears at tuq = 2.35 TeV, and the binding energy reaches the bootstrap point 
(E - 2m) /(2m) = -0.5 at m = 3.5 TeV. 

Despite the huge couplings required for binding, the corresponding physical mass is 
relatively small compared to that found in previous works |]-J|]. Fig.4 reports the binding 
energy versus m (physical mass) for A = 4 TeV. The onset of the bound-state solution is at 
m = 519 GeV, while the mass bootstrap point is reached at m — 386 GeV. We notice that 
the binding energy now increases with decreasing m. Moreover, according to Eq.flSBp, an 
even smaller mass is required for larger choices of the energy cut-off A. 

At the light of our study a Higgs-Higgs bound state would be conceivable for m ~ 
100 — 500 GeV provided that the coupling is very strong mo ~ 2 TeV. We stress the role 
played by mass renormalization in determining both, the shift of bound states towards 
higher coupling strengths, and the corresponding reduction of the physical mass required for 
bonding. Most of the previous calculations should be revised at the light of the present result 
in order to establish if mass renormalization has been correctly addressed. We just mention 
Rupp's Bethe-Salpeter approach where the chosen subtraction point gives m = m at 
any coupling. 
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FIGURES 



FIG. 1. Bubble-chain diagrams contributing to n^ oop in Eq. (|37|) 

FIG. 2. The physical Higgs mass m versus the bare mass mo (which fixes the coupling strength), 
for several choices of the energy cut-off A=2.8, 3.2, 3.6, 4.0, 4.4, 4.8, 5.2, 5.6, 6.0, 6.4, 6.8 and 7.2 
TeV. Data are also reported for the cross-over point from monotonic to non-monotonic behaviour 
occurring at A c =3.05 TeV. The dashed line represents the m = mo approximation which only holds 
in the weak-coupling regime mo < 0.3 TeV. 

FIG. 3. Binding energy E — 2m in units of 2m versus bare mass mo (coupling strength) for a 
cut-off A = 4 TeV. 

FIG. 4. Binding energy E — 2m in units of 2m versus physical Higgs mass m for a cut-off A = 4 
TeV. 
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